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The central-spin problem, in which an electron spin interacts with a nuclear spin bath, is a widely 
studied model of quantum decoherence (H- Dynamic nuclear polarization (DNP) occurs in central 
spin systems when electronic angular momentum is transferred to nuclear spins [3 and is exploited 
in spin-based quantum information processing for coherent electron and nuclear spin control Q]. 

However, the mechanisms limiting DNP remain only partially understood [^. Here, we show that 
spin-orbit coupling quenches DNP in a GaAs double quantum dot [^, even though spin-orbit cou¬ 
pling in GaAs is weak. Using Landau-Zener sweeps, we measure the dependence of the electron 
spin-flip probability on the strength and direction of in-plane magnetic field, allowing us to distin¬ 
guish effects of the spin-orbit and hyperfine interactions. To confirm our interpretation, we measure 
high-bandwidth correlations in the electron spin-flip probability and attain results consistent with 
a significant spin-orbit contribution. We observe that DNP is quenched when the spin-orbit compo¬ 
nent exceeds the hyperfine, in agreement with a theoretical model. Our results shed new light on the 
surprising competition between the spin-orbit and hyperfine interactions in central-spin systems. 


Dynamic nuclear polarization occurs in many con¬ 
densed matter systems, and is used for sensitivity en¬ 
hancement in nuclear magnetic resonance @ and for de¬ 
tecting and initializing solid-state nuclear spin qubits 
DNP also occurs in two-dimensional electron systems 
via the contact hyperfine interaction. In both self- 
assembled [SE3 and gate-defined quantum dots iQ- 
0, for example, DNP is exploited to create stabilized 
nuclear configurations for improved quantum information 
processing. Closed-loop feedback [l5| based on DNP, in 
particular, is a key-component in one- and two-qubit op¬ 
erations in singlet-triplet qubits i 113, [11. 

Despite the importance of DNP, it remains unclear 
what factors limit DNP efficiency in semiconductor spin 
qubits Q. In particular, the relationship between the 
spin-orbit and hyperfine interactions I9l-l21| has been 


overlooked in previous experimental studies of DNP in 
quantum dots. In this work we show that spin-orbit 
coupling competes with the hyperfine interaction and 
ultimately quenches DNP in a GaAs double quantum 
dot d, [H, even though the spin orbit length is much 
larger than the interdot spacing. We use Landau-Zener 
(LZ) sweeps to characterize the static and dynamic prop¬ 
erties of Asrit), the splitting between the singlet S and 
nis = I triplet T+, and the observed suppression of DNP 
agrees quantitatively with a new theoretical model. 


Figure [TTa) shows the double quantum dot used in 
this work [3, m. The detuning, e, between the dots 
determines the ground-state charge configuration, which 
is either (1,1) [one electron in each dot], or (0,2) [both 
electrons in the right dot] as shown in Fig.[T][b). To mea¬ 
sure Asrit), the electrons are initialized in |(0, 2)5'), e is 
swept through the S — T+ avoided crossing at e = esr, 
and the resulting spin state is measured [Fig.[2][a)]. In the 


absence of noise, slow sweeps cause transitions with near¬ 
unity probability. For large magnetic fields, however, we 
find maximum transition probabilities of approximately 
0.5. This reduction is a result of rapid fluctuations in 
the sweep rate arising from charge noise (see Supple¬ 
mentary Information). Even in the presence of noise, 
however, the average LZ probability (Plz(O) can be ap¬ 
proximated for fast sweeps as Supple¬ 

mentary Information). Here (■ • •) indicates an average 
over the hyperfine distribution and charge fiuctuations, 
and (3 = diEs — ET^)/dt is the sweep rate, with Es 
and Et_^ the energies of the S and r+ levels. To ac¬ 
curately measure asr = \J (l^sr(i)p), we therefore fit 
{Plz) vs (3~^ to a straight line for values of /3 such that 
0< (Plz) <0.1. [Fig.HKa)]. 

We first measure asr vs (j) a,t B = 0.5 T [Fig. [2][b)], 
where cj) is the angle between the magnetic field B and 
the z axis [Fig. [Ija)]. asr oscillates between its extreme 
values at 0° and 90° with a periodicity of 180°. Fixing 
(j) = 0° and varying B, we find that asr decreases weakly 
with with B, but when (p = 90°, asr increases steeply 
with B, reaching values greater than 10 times that for 
(/) = 0°, as shown in Fig. |2][c). 

We interpret these results by assuming that both 
the hyperfine and spin-orbit interactions contribute to 
Asrit) and by considering the charge configuration of 
the singlet state at esr [Figs.[T][b) and (c)]. The matrix 
element between S and r+ can be written as Asrit) = 
Anrit) + Aso- Anrit) = g*RB6B±it) is the hyperfine 
contribution, which arises from the difference in perpen¬ 
dicular (relative to B) hyperfine field, SB±it), between 
the two dots [i^. (In the following, we set g*R b = !■) 
Aupit), which is a complex number, couples [(1,1)5) 
to |(l,l)r+) when the two dots are symmetric. Aso 
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FIG. 1. Experimental setup, (a) Scanning electron micro¬ 
graph of the double quantum dot. A voltage difference be¬ 
tween the gates adjusts the detuning e between the potential 
wells, and a nearby quantum dot on the left senses the charge 
state of the double dot. The gate on the right couples the 
double dot to an adjacent double dot, which is unused in this 
work. The angle between B and the z axis is (f). (b) En¬ 
ergy level diagram showing the two-electron spin states and 
zoom-in of the S — T+ avoided crossing, (c) The hyperfine 
interaction couples |(1,1)S) and |(1, 1)7+) when the two dots 
are symmetric, regardless of the orientation of B, and the 
spin-orbit interaction couples |(0, 2)5') and |(1, l)r+) when B 
has a component perpendicular to flso = ^soz, the effective 
spin-orbit field experienced by the electrons during tunneling. 


is the spin-orbit contribution, which arises from an ef¬ 
fective magnetic field flso = flsoz experienced by the 
electron during tunneling [l9l |. Only the component of 
flso -L B causes an electron spin flip. Aso therefore 
couples 1(0,2)5) to |(1,1)T+) when cj) ^ 0°, and 
is proportional to the double-dot tunnel coupling |l9l |. 
which is 23.1 /reV here. At esr, the singlet state |5) is a 
hybridized mixture: |5) = cos0|(l, 1)5) -I- sin0|(0,2)5), 
where the singlet mixing angle B approaches 7r/2 as B 
increases (see Supplementary Information). Taking both 
B and (p into account, we write 0 

AsT(t) = Ansit) + Aso 

= SB±{t)cosB + flsosin(j)smB. (1) 

The data in Fig. [IDb) therefore reflect the dependence 
of Asrit) on p in equation ([T]). The data in Fig.[2Dc) re¬ 
flect the dependence of Asrit) on B. As B increases, 
B also increases, and |5) becomes more |(0,2)5)-like, 
causing Affrit) to decrease. When (/> = 0°, Aso = 0 
for all B, but when p = 90°, Aso = f^sosind, and 



(|) (degrees) B (T) 

FIG. 2. Measurements of asr- (a) Data for a series of LZ 
sweeps with varying rates, showing reduction in maximum 
probability due to charge noise. The horizontal axis is pro¬ 
portional to the sweep time. Upper inset: Data and linear 
fit for fast sweeps such that 0 < (Plz) <0.1. Lower Inset: 
In a LZ sweep, a |(0,2)5) state is prepared, and e is swept 
through esT (dashed line) with varying rates. Here h — 2nh is 
Planck’s constant, (b) gst vs p (dots) and simulation (solid 
line), (c) a ST vs B for p = 0° and p = 90° (dots) and fits to 
equation ([T]) (solid lines). Error bars are fit errors. 


asT increases with B. Fitting the data in Fig. [2Dc) al¬ 
lows a direct measurement of the spin-orbit and hyper- 
fine couplings (see Supplementary Information). We find 
yj{\5B'^{t)\) = 34± 1 neV and flso = 461 ± 10 neW cor¬ 
responding to a spin-orbit length Aso ~ 13 /rm IH 
good agreement with previous estimates in GaAs [22 

We further verify that Asrit) contains a signifi¬ 
cant spin-orbit contribution by measuring the dynam¬ 
ical properties of Przit)- A key difference between 
the spin-orbit and hyperfine components is that Aso 
is static, while Ah pit) varies in time because it arises 
from the transverse Overhauser field, which can be con¬ 
sidered a precessing nuclear polarization in the semi- 
classical limit 22|. To distinguish the components of 
Asrit) through their time-dependence, we develop a 
high-bandwidth technique to measure the power spec¬ 
trum of Ppzit)- 

Instead of measuring the two-electron spin state after 
a single sweep, e is swept twice through esr with a pause 
of length r between sweeps [Fig. |3Da)] (See Supplemen¬ 
tary Information). Assuming that Stiickelberg oscilla¬ 
tions rapidly dephase during r [isl . 0, and after sub¬ 
tracting a background and neglecting electron spin re- 
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FIG. 3. Correlations and power spectrum of PLzif). (a) 
Pulse sequence to measure Rpp(t) using two LZ sweeps, (b) 
Rpp{t) for (f> = ^° and B = 0.1 T. The data extend to r = 200 
/is, but for clarity are only shown to 75 /rs here, (c) Sp{ui) 
vs 4> obtained by Fourier-transforming Rpp(t). At </> = 0°, 
the differences between the nuclear Larmor frequencies are 
evident, but for \(l>\ > 0°, the absolute Larmor frequencies ap¬ 
pear, consistent with a spin-orbit contribution to ost- The 
reduction in frequency with 4> is likely due to the placement of 
the device slightly off-center in our magnet (see Supplemen¬ 
tary Information), (d) Line cuts of Sp{ uj) at <() = 0°, 25°, and 
80°. 


laxation, the time-averaged triplet return probability is 
proportional to Rpp{t) = {PLz{t)PLz{t + t )), the au¬ 
tocorrelation of the LZ probability [Fig. |3Kb)]. Taking a 
Fourier-transform therefore gives 5'p(w), the power spec¬ 
trum of PLz(t) [Figs. inKc) and[31[d)]. For Ppzit) <C 1, 

PLz{t) OC |A5T(t)|^ 

spectrum of |A 5 T(t)| 
us to measure the high-frequency components of Sp{uj), 
because the maximum bandwdith is not limited by the 
quantum dot readout time. 

Because it arises from the precessing transverse nuclear 
polarization, Appit) contains Fourier components at the 
Larmor frequencies of the ®®Ga, ^^Ga, and ^®As nuclei in 
the heterostructure, i.e., App(t) = , 

where a = 1, 2, or 3 indexes the three nuclear species, and 
the 9a are the phases of the nuclear fields. Without spin- 


so Sp{u!) OC 5'|Asr|2(w), the power 
^ This two-sweep technique allows 


orbit interaction, |A 5 T(t)|^ = \J2l=i 
contains only Fourier components at the differences of the 
nuclear Larmor frequencies. With a spin-orbit contribu¬ 
tion, however, |As 7 ’(t)p = IA 50 + App(t)|^ contains 
cross-terms like that give |A 5 T(t)p 

Fourier components at the absolute Larmor frequencies. 
A signature of the spin-orbit interaction would therefore 
be the presence of the absolute Larmor frequencies in 
Sp{uj) for [iil. 

Figure [SKb) shows Rpp{t) measured with i? = 0.1 T 
and (p = 0°. FigurelSKc) shows Sp{u!) for 0° < p < 90°. 
At (j) = 0°, only the differences between the Larmor fre¬ 
quencies are evident, but as p increases, the absolute nu¬ 
clear Larmor frequencies appear, as expected for a static 
spin-orbit contribution to Asrit)- These results, includ¬ 
ing the peak heights, which reflect isotopic abundances 
and relative hyperfine couplings, agree well with simula¬ 
tions (see Supplementary Information). 

Having established the importance of spin-orbit cou¬ 
pling at the S — T+ crossing, we next investigate how 
the spin-orbit interaction affects DNP. Previous research 
has shown that repeated LZ sweeps through esr increase 
both the average and differential nuclear longitudinal po¬ 
larization in double quantum dots Q. However, the rea¬ 
sons for left/right symmetry breaking, which is needed 
for differential DNP (dDNP), and the factors limiting 
DNP efficiency in general are only partially understood. 
Here, we measure dDNP precisely by measuring SB^, 
the differential Overhauser field, using rapid Hamiltonian 
learning strategies [i^ before and after 100 LZ sweeps to 
pump the nuclei with rates chosen such that {Plz) = 0.4 
(see Supplementary Information) [Fig. UKa)]. 

Figure HDb) plots the change in 5B^ per electron spin 
flip for B = 0.2 T and B = 0.8 T for varying p. In 
each case, the dDNP decreases with \(p\. Because the 
spin-orbit interaction allows electron spin flips without 
corresponding nuclear spin flops, dDNP is suppressed as 
j^soj = 1^50 sinsin01 increases with The reduc¬ 
tion in dDNP occurs more rapidly at 0.8 T because Aso 
is larger at 0.8 T than at 0.2 T. We gain further insight 
into this behavior by plotting the data against app / (Jst , 
where app = \J{\App{t)\‘^) [FigHDc)]. Plotted in this 
way, the two data sets show nearly identical behavior, 
suggesting that the size of the hyperfine interaction rela¬ 
tive to the total splitting primarily determines the DNP 
efficiency. 

Based on theoretical results and experimental data, to 
be presented elsewhere, we expect that the dDNP should 
be proportional to the total DNP, with a constant of 
proportionality that depends on B, but not /3 or </>. We 
therefore explain our measurements of dDNP using a the- 
oretcal model in which we have computed the average 
angular momentum {6m) transfered to the ensemble of 























4 


Probe Pump Probe 

6BzXl20 sweep X 100 6BzXl20 




(C) 



FIG. 4. DNP quenching by spin-orbit coupling, (a) Protocol 
to measure DNP. SBz is measured before and after 100 LZ 
sweeps by evolving the electrons around 5Bz- (b) dDNP vs 
(j> at fixed {Plz) = 0.4 for B = 0.8 T and B = 0.2 T and 
theoretical curves (solid lines). dDNP is suppressed for \(j)\ > 
0 because of spin-orbit coupling, (c) Data and theoretical 
curves for fixed {Plz) collapse when normalized and plotted 
vs (thfIost- Vertical error bars are statistical uncertainties, 
and horizontal error bars are fit errors. 


nuclear spins following a LZ sweep as: 

ISi \ 2 I PLzi^ST)\ 

{Sm)o,aHF^ lAg^l /- (2) 

where P'j^^{Ast) is the derivative of the LZ probabil¬ 
ity with respect to the magnitude of the splitting. (See 
Supplementary Information for more details.) Neglect¬ 
ing c harg e noise, we have the usual Landau-Zener for¬ 
mula 231 


Plz{Ast) = 1 — exp ( — 


27r|A5Th 


hp 


and equation @ reduces to 


27T 


{5m) oc (1 - Plz) ■ 


( 3 ) 


( 4 ) 


The data in Figs. SKb) and (c) can therefore be under¬ 
stood in light of equation dU because as the splitting ust 
increases with the sweep rate /3 was also increased to 
maintain a constant {Plz)- Because the hyperfine con¬ 
tribution uhf is independent of </>, {5m) therefore de¬ 
creases. The data collapse in Fig.SKc) can also be under¬ 
stood from equation ([4]), assuming a constant splitting 
and fixed probability. In this case, (3 oc jAsTp, as fol¬ 
lows from equation and hence {5m) oc a'^jp/\A st\^■ 
Measurements with fixed rate j3 also exhibit a similar 
suppression of dDNP (see Supplementary Information). 
In this case {Plz) increases with |(/)|, because of the in¬ 
creasing spin-orbit contribution to ost, and according to 
equation (HI, {5m) therefore decreases. 

The two theoretical curves in Figs.UKb) and (c) are cal¬ 
culated using equation O multiplied by fitting constants 
C, which are different for the two fields, and agree well 
with the data. As discussed in the Supplementary In¬ 
formation, we do not expect charge noise to modify the 
agreement between theory and data in Figs. SKb) and 
(c) beyond the experimental accuracy. Interestingly, the 
peak dDNP is less at B = 0.2 T than at B = 0.8 T, 
perhaps because the electron-nuclear coupling becomes 
increasingly asymmetric with respect to the center of the 
quantum dots at higher fields [28j. Finally, the peak 
dDNP value also approximately agrees with a simple cal¬ 
culation (see Supplementary Information) based on mea¬ 
sured properties of the double dot. 

In summary, we have used LZ sweeps to measure the 
S — Tp splitting in a GaAs double quantum dot. We find 
that the spin-orbit coupling dominates the hyperfine in¬ 
teraction and quenches DNP for a wide range of magnetic 
field strengths. A misalignment of B to Plso by only 5° 
at i? = I T can reduce the DNP rate by a factor of two, 
and DNP is completely suppressed for a misalignment 
of 15°. The techniques developed here are directly ap¬ 
plicable to other quantum systems such as InAs or InSb 
nanowires and SiGe quantum wells, where the spin-orbit 
and hyperfine interactions compete. On a fundamental 
level, our findings suggest avenues of exploration for im¬ 
proved S-T+ qubit operation [23| and underscore the 
importance of the spin-orbit interaction in the study of 
nuclear dark states and other mechanisms that 

limit DNP efficiency in central-spin systems. 
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1. MEASURING asr 

Here we describe the htting procedure to extract ctst- The experimentally measured 
quantity is the average triplet occupation probability (Pt), which we interpret as the aver¬ 
age Landau-Zener (LZ) probability (Plz), at the end of a sweep. Here (■ ■ ■) indicates an 
average over the hyperhne distribution and charge fluctuations for the same nominal swe^ 
parameters. We calibrate the rate fl = d{Es — Ej’j^)/dt using the spin-funnel technique [l| 
and assume a linear change in the S' — T+ splitting near the avoided crossing. 

^HF{t) varies in time because of the nuclear Larmor precession and statistical fluctua¬ 
tions in the magnitude of the nuclear polarizations. We argue that both types of hyperhne 
huctuations occur on time scales much longer than LZ transitions and can be treated as 
quasi-static. In typical experiments, the S — T+ splitting is swept through approximately 
5 GHz in less than 1 ps. For splittings of order 10 MHz, the total time spent near the 
avoided crossing is less than 10 ns, which is much faster than the nuclear Larmor period at 
1 T, roughly 100 ns. Furthermore, during 1 /rs, the nuclear polarization diffuses by approx¬ 
imately 7 kHz j^, which is 3 orders of magnitude smaller than anp. We therefore assume 
that the splitting is constant during a single sweep. Numerical simulations discussed below 
also support the hypothesis that nuclear Larmor precession does not signihcantly affect (Pt) 
for the sweep rates used here [Fig. IST] . 

In the absence of hyperhne or charge huctuations, the probability for a transition is given 



2 


by the LZ formula: Pizif) = 1 — exp(—27r|A5'r(f)p/(h/9)) [31]. Neglecting high-frequency 
charge noise, the exact form of the LZ probability averaged over the hyperhne distribution 
can be computed. Let the total splitting be ^st = ^hf + ^so- We take A 50 to be the 
constant, real spin-orbit part and /S.hf the complex hyperhne contribution. Assuming that 
the real and imaginary parts of ^hf {u and v, respectively) are Gaussian-distributed around 
zero such that the root-mean-square hyperhne splitting is anF, the probability distribution 
for the splitting to have magnitude A = jA^'p] is 


p(A) = 


du 


dv 


TTCr 


HF J-00 


2 I 2 

_ u ~h'^ / \ 

e '^HF 6 yA- (A50 + j 


2A - 

e 




a 


^hf Iq{2AAso/(Thf)^ 


(51) 

(52) 


HF 


where Jq is the zeroth-order modihed Bessel function of the hrst kind. Note that when Aso = 
0, equation flS^ reduces to the familiar distribution p(A) = ^"^hf j^j. Integrating 

the LZ probability over this distribution yields the average LZ probability {Plz)- 

27rA2' 


{Plz) = J^ dA(^l-exp 


hfd 


p{A) 


= 1 — Q exp 




Q 


(53) 

(54) 


with 


Q = 




hi3 


fl. 


(S5) 


Note also that to leading order in 


Note that this result agrees with another derivation 
/5"^ {Plz) ~ 27r (A|o + ajjp) /h/d. 

The average triplet return probability (Pt) may be modihed due to ehects of charge 
noise on the dehning gates or in the two-dimensional electron gas itself. High-frequency 
charge noise in double quantum dots has recently been identihed as a major source of 
decoherence [61]. In the current setting, corrections to (Pt) should occur, because charge 
huctuations lead to time-dependent variations in S' — T+ detuning Es — i?T+, on top of the 
linear time-dependence due to the prescribed sweep rate fd. Additionally, charge huctuations 
can add noise to the oh-diagonal coupling Asrit) = 6B±(t) cos 9 + Gsosin^sin^, because 


the singlet mixing angle 6 = tan 


2t 


depends on e. (Here t = 23.1 peV is the 
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B=0.1 T, (})=90° 



B=0.3 T, (t)=90° B=0.5 T, (l)=90° B=0.7 T, (|)=90° B=0.9 T, (|)=90° 




Data Simulation with Simulation with Analytic 

charge noise, hyperfine averaging, 

hyperfine averaging, Larmor precession 

Larmor precession 


FIG. SI. Comparison of LZ data and simulations. Each panel shows data and simulations for a 
different magnetic field strength and orientation. Red curves are experimental data for a series of 
LZ sweeps with varying rates. Blue curves are simulated data including charge noise, hyperfine 
averaging, and nuclear Larmor precession for the calculated value of the splitting corresponding 
to the red curves. Green curves are simulated data with hyperfine averaging and Larmor pre¬ 
cession for the same value of the splitting as the blue curves. Black curves are calculated via 
equation (ISill using the same value of the splitting. In all panels, the y axis is {Plz), and the x 
axis is h/f3 (/is/GHz). Here h = 2Trh is Planck’s constant. 


double-dot tunnel coupling.) As discussed below, however, the noise in Ast should have 
much less effect than the detuning noise for the magnetic helds studied here. 

We observe that for high magnetic helds and slow sweeps, the maximum LZ probability 
falls to 0.5 as shown in Fig. EH It was previously noted that strong detuning noise can 
have such an effect j^. To conhrm that charge noise causes the probability reduction, we 
have performed Monte Carlo simulations of the Schrodinger equation for symmetric double 
dots undergoing LZ sweeps, including the effects of wide-band charge noise, nuclear Larmor 
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precession, and averaging over the hyperfine distribution. The results of the simulations 
and experimental data are shown in Fig. [STl We generate random charge noise with power 
1-14 { iH^V ^ 


spectrum 14 x 10“^^ ^ for / < 1 GHz, and 0 otherwise. We generate the Fourier 

transform of the charge noise time record by picking the amplitude corresponding to the 
chosen power spectrum and a random phase for each frequency / in the desired range. We 
then perform an inverse Fourier transform to obtain the charge noise time record. The 
spectrum we chose corresponds to a noise amplitude of 3 nV/\/Hz at / = 1 MHz, which is 
approximately the measured level of charge noise in the double dot used here. Note that we 


have extrapo 
MHz in ref. 


ated the frequency dependence that was previously measured to / = 1 
6| up to / = 1 GHz in these simulations. However, one expects the results 
to be most sensitive to noise in the range of 10-100 MHz, corresponding to the size of the 
splitting. The e-dependent Hamiltonian used in these simulations was 


H(e) = 




, 6B1 (t) cos 6 + Qso sin (j) sin 6 


SB±(t) cos 9 + flso sin 0 sin 9 


(S6) 


in the {|T+), 15)} basis. Linear e sweeps through the 5 — T+ crossing esr = were used 
in the simulation to replicate the actual experiments. For each strength and orientation of 
the magnetic held, 9 was calculated at esT using the measured tunnel coupling, and the htted 
values of the spin-orbit and hyperhne couplings from the main text were used to compute 
the splitting. We assumed a lever arm of 10 to convert the voltage noise on the quantum 
dot gates to e noise. 

The simulated LZ curves with charge noise agree well with the data as shown in Fig. [STl 
The same simulations including averaging over the hyperhne distribution and nuclear Larmor 
precession, but without charge noise, show very little reduction in probability compared with 
the analytic result, equation fIMl) . supporting the hypothesis that charge noise is responsible 
for most of the observed probability reduction. A key feature in these experiments is the 
decreasing maximum probability with increasing magnetic held. We can understand that 
this trend occurs because the ehect of charge noise on the Landau Zener probability is 
controlled by the huctuation in the energy splitting 6E{e) produced by a given huctuation 
in the detuning e, which is proportional to Since E{e) = ^ — B + |\/e^ 4- 4F, the 

magnitude of increases sharply with increasing magnetic held. 
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FIG. S2. Fitting error. We compute the fitting error by simulating {Plz) for the calculated splitting 
at each of the magnetic field configurations in the presence of charge noise. The simulated {Plz) vs 
is fitted to a straight line for 0 < {Plz) < 0.1, and the fitted value of the splitting is subtracted 
from the value chosen for the simulation. The difference is then divided by the simulated value of 
the splitting. Error bars are fit errors. 

Even in the presence of noise, however, the average LZ probability in the limit of fast 
sweeps is still 27r|A5'r(t)p/h/3, which is identical to the leading order behavior of the usnal 
LZ formula, as shown in section 3.1 of ref. [^. Replacing the LZ formula in equation flS4p by 
its leading order behavior, and performing the integration over the quasi-static distribution 
gives {Plz) ~ 27r (A|q -|- a'jjp) /h(3. Such a result can be understood because the effect of 
detuning noise is reduced on short time scales. Figure ISTl demonstrates this idea because the 
analytic curves deviate significantly from the data for {Plz) ^ 0.2, but for 0 < {Plz) < 0.1, 
the analytic results agrees well with the data. Based on additional simulations, we estimate 
the systematic error in the deduced value of asr as obtained by fitting measured values of 
{Plz) to a straight line for values of (3 such that 0 < {Plz) < 0.1 to be small for most of 
the experimental conditions as shown in Fig. [S^ 

We note that the coupling Asrit) = 6B±{t) cos 9+flso sin 0sin 9 depends on e through the 
singlet mixing angle 9. This dependence means that during a LZ sweep, the coupling Asrit) 
varies both due to the linear e sweep as well as charge noise. We estimate that > 40 ^^^: 
for the fields studied here. We therefore expect detuning fluctuations to be the dominant 
noise source. Furthermore, when \E{e)\ < asr, o'st changes by only a few percent during 
the sweep and is likely not a significant source of error in the measurement of Asrit)- 
Additionally, we note that the simulations in Fig. [STl which include e-depending coupling, 
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demonstrate that the htting procedure described above allows an accurate measurement 
of (Tst- Finally, we have also performed additional simulations, taking into account the 
measured values of E{e), which deviate slightly from the values predicted by assuming a 
constant tunnel coupling, and we observe no signihcant change in our results. 


2. DIRECTION OF ft 


so 


The double quantum dot axis is aligned within ^ 5° of either the [110] or [110] axes 
of the crystal, but we do not know which. In the later case, both the Rashba and Dres- 
selhaus spin-orbit helds are aligned with the z axis, and their magnitudes add 8|. In the 
former case, the Rashba and Dresselhaus contributions are also aligned with the z axis, but 
their magnitudes subtract. The techniques used here could be employed to distinguish the 
Rashba and Dresselhaus spin-orbit contributions by measuring asT with double quantum 
dots fabricated on different directions with respect to the crystal axes. 


3. FITTING asT VS B AND 


We £t the data in Fig. 2(c) in the main text to a function of the form asT = 


V^so sin^ 6 sin^ cj) + ajjp cos^ 6, with Aso and anp as £t parameters. The singlet mixing 
angle 9 is computed by assuming that the (1,1) and (0,2) singlet branches are a two-level 
system with constant tunnel coupling, as discussed above. 

Aso is held at 0 when htting data for 0 = 0° to determine the hyperhne coupling. We 
also exclude data points for B < 0.2 T in the Ft, as the hyperhne contribution appears to 
decrease at very low helds. We determine the spin orbit length using equation (28) of Ref. 
8], where the spin-orbit held is computed as ilso = y ^Xso ’ ^dqd ~ 200 nm is the 

interdot spacing, and Xso is the spin-orbit length. The simulation in Fig. 2(b) in the main 
text is generated using the same equation with the htted values of the Aso and app- 


4. MEASURING Rpp{t) 

Here we derive the triplet return probability after two consecutive LZ sweeps with a 
pause of length r in between. In experiments, both sweeps were in the same direction, and 
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e was held in the (0, 2) region between sweeps, as shown in Fig. 3(a) in the main text. 
Suppose the hrst LZ sweep takes place at time t with probability PLzif)- The probability 
for the two electrons to be in the T+ state is PLzif), while the probability to be in the 
S state is 1 — Pizif)- Then, the detuning is quickly swept into the (0,2) region. Here, 
electron spin dephasing occurs rapidly, and there is very little T+ occupation in thermal 
equilibrium because the S and T+ states are widely separated in energy. Thus, after a wait 
of length r, but before the second sweep, the triplet population is PLz{t)e~'^^'^^ , and the 
singlet population is 1 — PLz{'t)e~'^^'^^, where Ti is the electron relaxation time. After the 
second sweep, the triplet occupation probability is 

Prit + r) = (1 - PLz{t)e-^l^^) PLz{t + r) + PLz{t)e-^^^^ (1 - PLz{t + r)) (S7) 

= -2PLz{t)PLz{t + r)e-"/^i + PLz{t + r) + ■ (S8) 


The second and third terms in equation flS8D vary slowly with r. These terms are found by 
htting the measured triplet probability to an exponential with an offset and are subtracted. 
When T\ 3> r, relaxation can be neglected, and the predicted time-averaged signal is (Pt(^ + 
r)) cx i?pp(r), where Ppp(r) = (Pp^(t)Pp 2 (t-|-r)), the autocorrelation of the LZ probability. 
When 0 = 0°, Ti 3> r^ax = 200 /rs, where T^ax is the largest value of r measured. The 
shortest relaxation time Ti ~ 100 /is in these exj^riments time occurs when 0 = 90°, which 
is consistent with spin-orbit-induced relaxation 9(]. 

The effect of Ti relaxation is to multiply the measured correlation by an exponentially- 
decaying window, which reduces the spectral resolution of the Fourier transform, but does 
not shift the frequency of the observed peaks. We expect statistical fluctuations in the 
amplitude of the hyperfine field to affect the spectrum in a similar way, although we expect 
this effect to be less than that of electron relaxation. The raw data, [Fig. 3(b) in the 
main text] consisting of 667 points (each a result of two sweeps with a 40 % chance of a 
LZ transition) spaced by 300 ns, were zero-padded to a size of 1691 points to smooth the 
spectrum, and a Gaussian window with time constant 150 /is was applied to reduce the 
effects of noise and ringing from zero-padding before Fourier transforming. 

The magnetic resonance frequencies in Fig. 3(c) appear to decrease with 0. The inhomo¬ 
geneity of the x-coil in our vector magnet is 1.6 % at 0.6 cm offset from the center. Thus, 
the field could easily be reduced by more than 3 % for a misplacement of the sample by 1 
cm from the magnet center. We have simulated the data in Fig. 3(c) in the main text based 
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FIG. S3. Simulations of Sp{u}). (a) Experimental data, (b) Theoretical simulation taking into 
account known sweep rates, nuclear magnetic resonance frequencies, hyperfine couplings, and a 
4.4% reduction in field in the x direction. The expected frequencies at S = 0.1 T are f69Qa = 
1.0248 MHz, f7iGa = 1-302 MHz, and fr^As = 0.7315 MHz. 


on the measured hyperfine and spin-orbit couplings and the known sweep rates. Assuming a 
4.4 % reduction in the held from the x-coil, we obtain good agreement between theory and 
experiment [Figs. [S3f a) and (b)]. 

We argued in the main text that only the difference frequencies should appear in the spec¬ 
trum Sp{u) without spin-orbit coupling by considering the time-dependence of |As'r(f)p and 
because Sp{u) oc when Pizif) "C 1. Since Pizif) contains only even powers of 

|A 5 r(f)|, Sp{uj) can generally be expressed in terms of differences of the resonance frequen¬ 
cies, but will not contain the absolute frequencies in the absence of spin-orbit coupling, 
regardless of the value of Ppzii)- 
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5. DERIVATION OF NUCLEAR POLARIZATION CHANGE {5m) 


Here we derive equations 2 and 4 in the main text. Let ^st = ^so ^ ^hf where ^so 
is real and 

A/+, (S9) 

j 

where is the raising operator for the nuclear spin, and the Xj are individual coupling 
constants. We assume that there are many nuclear spins, so that each coupling constant is 
small. Also, 

^ (I = ?/(/ + 1) = I E (SIO) 

j i 

where / = | is the spin of the nuclei, and the angular brackets refer to an average over the 
distribution of nuclear spins. 

We pick one of the nuclear spins, j, and we wish to compute {5mj), the mean value of the 
change in /| after one sweep. Let P^zi^sr) be the probability of an S' — T+ transition for 
a hxed value of Ahf- Clearly, Plz depends on |A 5 r|- We calculate Srrij as follows. Write 

AsT = a + bP^\ (Sll) 


where a includes the contributions of spin orbit and of all nuclei other than the nucleus j, 
and the second term rraresents the contribution (of order Xj) from nucleus j. According to 
equation (31) of Ref. the value of 5mj for this conhguration should be given by 


drrij = — 
^ 27r 


dOjPiziAsT) 


dip 

dFj' 


(S12) 


where p = arctan(Im(A 5 T)/Re(A 5 T)) specihes the orientation of Ast in the complex plane. 
Without loss of generality, we may suppose that a is real. Then we have, ignoring terms 
that are higher order in 6/a, 

rJfn h 

(513) 

(514) 


dp b 
— = -cosOj 
dOj a 


Plz{Ast) — Plz{o) + bP'i^zi^) cos6j 
drrij = ^Pizia), 


(S15) 


where P'hzi^) is the derivative of Plz{o) with respect to a. Averaging over nuclear conhgu- 
rations, we obtain 


{drrij) = (6^) 


P'lz^^) 

2a 


(S16) 
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with (6^) = (5/2) A^. In the case of no charge noise, we have 


Plz{^st) = 1 - exp 


2vr|A5'rp 

h(3 


so 


Plz{^) 

2a 




and 




LZ 


»)• 


(S17) 


(S18) 


(S19) 


Finally, we sum over all nuclear spins and make the replacement a 

2tt 

- Plz{^st))- 


[Asrl, obtaining 


(S20) 


The collapse demonstrated in Fig. 4(c) in the main text can be understood from equa¬ 
tion flS20p . assuming constant ^st and hxed probability. In this case, (3 oc lA^rp from 
equation (1S17D . and hence {5m) oc Ag-pp. 

In the case of a fixed splitting, equation flS20p reduces to 

27r 


(<5m) = exp 


27r|A5rP 


(S21) 


In equation flS2ip . {5m) —)■ 0 for both /3 —)■ 0 and (3 ^ oo. In practice however, experiments 
necessarily average over the hyperhne distribution. Thus, using equation flS4|) with Aso = 0 
to compute (1 — Plz{.5S.st)) i we have 

2f 


= -^^hfQ 


27rcrl 


hi3 




(522) 

(523) 


hp 


According to equation flS23p . in the limit of slow sleeps, where /5 —)■ 0, {5m) —)■ 1, and in the 
limit of fast sweeps, where (3 ^ oo, {5m) —)■ 0, as expected. 

The theory curves in Figs. 4(c) and (d) in the main text were generated by computing 
equation flS20p . For each held angle (j), the parameters 9, Aso, and anF were calculated 
using the htted values of the spin-orbit and hyperhne couplings as well as the measured 
tunnel coupling. Equation flSip was then solved using the calculated parameters to hnd 
the rate [3 such that {Plz) = 0.4. In order to compare with data on the dDNP rate, the 
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FIG. S4. DNP quenching with fixed sweep rate, (a) dDNP ys (f) at B = 0.2 T and B = 0.8 T. For 
each held, the sweep rate /3 was chosen to give {Plz) = 0.4 at 4> = 0° and then was held constant 
for (j) 7 ^ 0°. (b) As \(j)\ increases, ctst increases. As a result, {Plz) also increases and DNP is 
suppressed, according to equation (|S20I) . Error bars are statistical uncertainties. Lines between 
points serve as a guide to the eye. 


theoretical curves for {5m) were multiplied by fitting constants G, which are different for the 
two curves. As explained in the main text, and further discussed below, we expect the ratio 
between the dDNP rate and {5m) to depend on the magnetic held but to be independent of 
the sweep rate. 


Data taken at hxed sweep rate (3 also show a suppression of DNP, as shown in Fig. [Mf ai. 
In this case, {Plz) increases with 101 because of spin-orbit coupling [Fig. ISlf bi]. and {Plz) 
therefore increases, causing {5m) to decrease, according to equation flS20|) . 


To address the effect of charge noise on dDNP, we recompute equation flS20p in the limit of 


strong noise using the results of Ref. j^, making the replacement P{a) = ^ (1 — exp 


hi3 


for Plz{o) both in the derivation leading to equation (IS20p and in equation flS4l) for the 


computation of f3. The expected dDNP in the presence of strong noise is shown in Fig. [S5l 
and it does not signihcantly deviate from the case without noise, at least at the level of the 
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FIG. S5. The effect of charge noise on dDNP. The data and solid lines are the same as in Fig. 4 
in the main text, and the dashed lines are the theoretical estimates for dDNP in the presence of 
charge noise. The dashed and solid lines are normalized to the same values at (f) = 0°. Error bars 
are statistical uncertainties. 

experimental accuracy. 


6. MEASURING 5B, 


We measure SB^ by first initializing the double dot in the |(0, 2)^) state and then separat¬ 
ing the electrons by rapidly changing e to a large negative value [1^. When the electrons are 
separated, the exchange energy is negligible, and the magnetic held gradient 6Bz drives oscil¬ 
lations between \S) and \Tq). In our experiments, we measure the two-electron spin state for 
120 linearly increasing values of the separation time. The resulting single-shot measurement 
record is thresholded, zero padded, and Fourier transformed. The frequency corresponding 
to the peak in the resulting Fourier transform is chosen as the value of SB^. This technique 
is related to a previously described rapid Hamiltonian estimation technique [^. 


7. EXPECTED DNP RATE 

In this section we give a simple calculation to explain the value of the peak (0 = 0°) 
dDNP rate, as shown in Fig. 4 of the main text. Additional measurements were carried 
out to measure the pumping rate of the sum hyperhne held, {Br + Bi)/2, where B^ and Bi 
denote the longitudinal hyperhne helds in the right and left dots. This rate was determined 
by measuring the location of esT before and after a series of LZ sweeps to polarize the nuclei 
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aX B = 0.2 T. We observe that the sum held is pumped roughly twice as efficiently as the 
difference field, = Bj. — Bi. Setting {Br + Bi)/2 = 2{Br — Bi), where indicates 
the pumping rate of the left(right) dot, we have Bi = (3/5)5^, meaning that the left dot is 
pumped 3/5 as often as the right dot. Under these conditions, the average gradient builds 
up at a rate (per electron spin flip) of {Br — Bi)/{Br + Bi) that is only 1/4 the rate that 
would occur if nuclear spin hips occurred in only one dot. 


To determine the expected change in SB^, we require the approximate number of spins 
overlapped by the electronic wave function in the double dot. We have measured the inho¬ 
mogeneous dephasing time of electronic oscillations around SB^ and find = 18 ns 10|. 
This dephasing time corresponds to a rms value of the gradient (Tsb^ = \/{\5Bz\^) = 
h/ (|5f*|/rB\/27rT2*) = 2 mT, where h is Planck’s constant. The total number of spins over¬ 
lapped by the wavefunction is iV = (hi/cx^s^)^ ~ 3 x 10®, where hi = 4.0 T 11]. If all nuclear 
spins were fully polarized, then the dots would experience a hyperfine field of ho = 5.3 T [n|. 
and if the nuclear spins in the two dots were fully polarized in opposite directions, the gra¬ 
dient would be 2ho. Therefore, the expected change in the gradient per electron spin flip, 
corresponding to a change in nuclear angular momentum of h, is ^ x ^ 2 j(Af/ 2 )° “ kHz, 
where I = 3/2 is the nuclear spin. The average dDNP under actual conditions is 1/4 of this 
value, or 3 kHz, in reasonable agreement with our observations. In addition, we note the 
reasonable agreement between the measured value of =2 mT and the root-mean-square 
hyperfine gap {\6B±{t)\‘^) Ri 34 neV/ = 1-5 mT. 
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